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Asymmetric Wolf-Rayet winds: implications for GRB 
afterglows. 
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ABSTRACT 

Recent observations of Wolf-Rayet (WR) binaries WR151 and WR155 infer that their 
stellar winds are asymmetric. We show that such asymmetries can alter the stellar- 
wind bubble structure, bringing the wind-termination shock closer to the WR star. 
If the wind asymmetry is caused by rotation, the wind density and distance to the 
wind-termination shock are both decreased along the rotation axis by a factor of a few 
for the observed equator-to-pole wind density ratio of WR151. If this asymmetry lasts 
until core-collapse the time taken to reach the wind-termination shock by supernova 
ejecta or a gamma-ray burst jet is reduced. This leads to a distorted structure of the 
supernova ejecta and makes it more likely a constant density environment is inferred 
from gamma-ray burst afterglow observations. 

Key words: stars: Wolf-Rayet - circumstellar matter - stars: individual: WR151 - 
stars: winds, outflows - gamma-rays: bursts 



1 INTRODUCTION 

Wolf-Rayet (WR) stars are the most massive stars in the 
final stages of evolution and have lost (or in the process of 
losing the last remnants of) their hydrogen envelopes. Ob- 
servationally they are very luminous, with broad emission 
lines indicating a fast and dense wind. WR stars are the 
preferred progenitors of type Ibc supernovae (SNe) due to 
the lack of a hydrogen envelope. They have also become 
the leading candidates to be the progenitors of Gamma-ray 
bursts (GRB). A GRB is thought to occur at core-collapse 
if a black hole forms and the core material surrounding the 
black hole possesses enough angular momentum to reside in 
an accretion disk. This disk can then feed material to the 
black hole resulting in a highly relativistic jet. If the progen- 
itor is compact the jet will emerge at the surface to produce 
the prompt emission and later the afterglow. The supporting 
evidence is that GRBs occur in star forming regions of their 
host galaxies, some GRBs have associated type Ibc SNe and 
some afterglow lightcurves are consistent with the jet propa- 
gating in a free-wind enviro nment expected around m assive 
stars (for a recent review see IWooslev fc Bloom! (|2006j). The 
evidence for stellar progenitors from absorption lines in af- 
terglow spectra h as weaken due to the recent findings of 
IChen et al.l (|2006l ). 

However when the circumburst environment (CBE) 
density structure is inferred from afterglow observations 



E-mail: j.eldridge@qub.ac.uk 



some GRB are found to occur in constant density en- 
vironment s rather than a freely expanding wind where 
o oc 1/r 2 llChevalier, Li fc FranssorJ 12004 |Panaitescull2005l : 
Ivan Marie et al. 20061 ). Stellar- wind bubbles do have a re- 
gion where the wind density structure is approximately a 
constant density. The expanding wind goes through a wind- 
termination shock (a hydraulic jump) when the ram pres- 
sure of the fr ee- wind drops below the pressure of the stalled- 
wind region (|Castor. McCrav fc Weaverll 19751 ). This stalled- 
wind region has a density profile that is roughly constant. 
Therefore if the free-wind region is small, the GRB jet can 
traverse it before the afterglow emission begins and a con- 
stant density environment would be inferred. Simulations 
show it is very difficult to shrink the free-wind region to a 
small enou gh radius for it to remain unobserved in the GRB 
afterglow (|van Marie et al.ll2006l : lEldridge et al.1 12006). The 
required radius is roughly 10 1 8 cm, a fraction of a parsec 
l|Chevalier. Li fc Fransson1l2004l ) . There are a number of pos- 
sible physical proc esses that might s hrink the free- wind re- 
gion as outlined by Ivan Marie et al. I l|2006l) but factors must 
be combined to produce a free- wind region with the required 
radius. 

One possible process that may have an effect is asym- 
metry in WR star winds. The cause of asymmetric winds 
would typically be the result of stellar rotation or du- 
plicity, however, as we discuss below there is some uncer- 
tainty in predicting how rotation affects the WR winds. 
Wind asymmetry can be estimated by performing po- 
larimet ry on stars as asymmetry indu ces strong polari- 
sation (|Brown. McLean fc Emslig Il978l ). Polarimetry sur- 
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veys of WR stars have been performed and for most 
WR stars the polarisation is low with only 20 percent 
of the stars showing any noticea ble intrinsic polarisation 
l|Harries. Hillier fc Howarthl [l998) . The observations indi- 
cate a mass-loss asymmetry between the highest and lowest 
wind density to be between a factor of 2 to 3. There is one 
important case of stronger polarisation from the close WR 
binary WR151 jVillar-Sbafn et alj|2006h . In this system the 
light is polarised by 4 percent and modelling indicates that 
the eqnator-to-pole density r atio (a) is 5. 

IVillar-Sbaffi et ah! (|2006l ) suggest that the difference is 
due to rotation redistributing the mass-loss from the pole 
to the eq uatorial region as in the wind - comp ressed zone 
model of llgnace. Cassinelli fc Biorkmanl (|199 6) . If this is 
the case then there is a direct implication for the CBE of 
GRB progenitors which are thought to be rapidly rotat- 
ing WR stars but the effect of stellar rotation on stellar 
winds is somewhat unc ertain. The effect on B star winds 
has been well studied. iBiorkman fc Cassine H (|l993h and 
llgnace. Cassinelli fc Biorkman 1 19961 ) first suggested that 
inertial effects due to rotation reduce the polar wind den- 
sity and increase the equatorial wind density. A more rie 
orous approach from Owocki. Cranmer fc Blondin ( 1994 



lOwocki. Cranmer fc Gavlevl (|l996l ) and iPetrenz fc Pulsl 
(2000) suggest the converse is true when the distortion of 
the stellar surface by rotation and the radiative driving force 

are taken into account; 

However only llgnace. Cassinelli fc Biorkmanl (| 19961 ) 
have considered WR winds. WR winds are different from 
winds of other hot stars as they are optically thick and hy- 
drodynamic models o f the winds are quite comp lex and com- 
putational intensive l|Grfener fc Hamannll2005r ). The results 
for B star winds cannot be directly applied to WR winds be- 
cause although they are both line driven the optical thick- 
ness may cause important differences and possibly result 
in the inertial terms domin ating and causing the winds to 
be dis torted as predicted by llgnace. Cassinelli fc Biorkmanl 
<| 19961 ). 

In this letter we use the observed asymmetry of 
WR151's wind to produce a model wind bubble for this sys- 
tem and study the effect which the asymmetry has on the 
bubble structure. However we do not consider what might 
cause the asymmetry. This avoids the uncertainty in the dif- 
ferent models for wind distortion by rotation. We then per- 
form a small parameter survey varying the initial interstel- 
lar medium (ISM) density and the size of the wind-density 
asymmetry and estimate the time it would take a GRB jet 
to reach the stalled-wind region. 



2 THE CURRENT AND FUTURE STELLAR 
WIND BUBBLE OF WR151. 

First to model the stellar- wind bubble we must select a stel- 
lar model to use. In WR151 the current mass of the WR star 
is estimated to be 23Mq. To obtain the correct mass from 
our single star model we adopt a solar metallicity (Z = 0.02) 
star with an initial mass of 6OM0 that has a long WN phase. 
We therefore assume that WR151 was initially a wide binary 
where the 6OM0 star became a red supergiant, filling its 
Roche Lobe and growing to a radius greater than the binary 
orbit so that common envelope evolution occurs. During this 



the orbit shrank to its current size, the hydrogen envelope 
was lost and the helium core became a WN star. For the 
secondary, we assume that the primary's hydrogen envelope 
was lost rapidly and the secondary accreted very little mate- 
rial from the primary. Therefore we assume that the initial 
mass of the secondary was similar to its current mass of 
3OM0. When we compare the wind strengths we find the 
primary wind dominates so we ignore the secondary's winds 
when modelling the wind bubble. 

Our stellar-win d bubble mode l is ma de using the same 
code as described in lEldridge et al.l l|2006l ) . We use the mass- 
loss rates and wind velocities from the solar metallicity 
6OM0model with an initial ISM density of no = 10cm -3 . 
We turn on the wind asymmetry when the star has lost its 
hydrogen envelope and assume that M(8) oc 1 + /3(sin#) 2 
where beta is j3 = a — 1. We normalise the mass-loss rate 
so that the total mass-loss from the star is unchanged from 
the symmetric case. Also we assume that the wind speed 
is unchanged. We select the stellar-wind bubble model with 
an age of 3.8 Myrs as the mass-loss rate and wind velocity 
are broadly in agreement with the observed parameters of 
WR151. 

Figure [T] shows a greyscale plot of the wind bubble den- 
sity structure with and without the asymmetric mass-loss. 
The dense shell in the stalled-wind region is the RSG wind 
that has been swept up by the faster WR wind. The differ- 
ence between the symmetric and asymmetric cases is clear 
with the distance to the wind termination shock in the po- 
lar direction reduced by nearly a factor of two. This is due 
to the reduction in wind density in the polar direction and 
therefore a reduction in the ram pressure compared to the 
stalled-wind ambient pressure which remains constant. 

Figure [2] compares density profiles through this envi- 
ronment. In this plot it is easier to compare the change 
in the position of the wind termination shock in the polar 
and equatorial directions. It is important to note that while 
the polar distance is decreased the equatorial distance to 
the wind termination shock is increased because the slightly 
higher density of the wind leads to a greater ram pressure. 
Also in these plots the asymmetry we have introduced has 
only just been switched on within the simulation and there- 
fore the stalled wind region is similar between the two sim- 
ulations. 

If we continue the simulations for another 0.34 Myrs 
up until the end of the stellar models (about a year before 
core-collapse) and we assume the asymmetry remains the 
distortion increases as the stalled-wind region is also altered. 
This is important as the supernova or gamma-ray burst that 
may occur at that time will propagate through this environ- 
ment. The resulting wind-bubble structures are shown in 
Figures [3] and [4] Both with and without the asymmetry the 
stalled-wind region has some denser swirls. These are due to 
insta bilities in the wind (|Garcia-Segura. Langer fc Mac-Low! 
Il99rj ) th at have grown to produ ce great density variations as 
found by lEldridge et al l (2006). We find with the asymmet- 
ric wind acting for a long time the magnitude and number 
of these swirls are greater. Overall we find the distance to 
the wind-termination shock, Rsw is decreased in the polar 
direction by a factor of approximately 3. 

If a SN was to occur in such a medium the ejecta 
would reach the wind-termination shock in the polar direc- 
tion sooner than along the equatorial region and therefore 
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Figure 1. Greyscale plots of the stellar-wind bubble structure within 12.5 parsecs of the WR star. The WR star had an initial mass of 
60Mq and is pictured here at an age of 3.8Myrs, the initial ISM density was 10cm -3 . The left panel shows the case for a spherically 
symmetric wind while the right panel shows the case when the equator-to-pole wind density ratio, a = 5. The polar direction is along 
the x-axis in both figures. 



decelerate more quickly by accreting more material. The SN 
remnant would be distorted to an oblate shape rather than 
being spherically symmetric. If a GRB was to occur in such 
a structure the afterglow jet would reach the stalled- wind re- 
gion in a shorter time and therefore the afterglow would ap- 
pear to be traversing a constant density environment. This 
is in comparison to the symmetric wind bubble case where 
a freely expanding wind would be inferred. 

There is an other reason why the -Rsw is more un- 
even in Figure EO Eldridgc et al. (2006) found that the po- 
sition of the wind-termination shock is not alway directly 
re lated to the initial ISM density a s the analytic relation 
of ICastor. McCrav fc Weaverl (|l975l ). This is because if the 
wind speed and wind density of the source star are vary- 
ing on the same timescale as the time it takes for the wind 
to reach the wind termination shock then the boundary will 
move. As that time now varies between the pole and equator 
so too does the magnitude of this timescale effect. 



3 A PARAMETER SEARCH. 

We have performed a parameter search varying the wind 
asymmetry, a, and the initial ISM density, no, around the 
star to investigate the resulting change in the polar dis- 
tance to the wind-termination shock. In table [1] we show 
the distance in parsecs to the wind termination shock and 
the wind density along the polar direction. Clearly the effect 
of increasing the wind asymmetry is that both Rsw and A* 
decrease. The effect of increasing the initial ISM density de- 
creases Rsw, as expected. There is one exception when at 
the highest ISM density, Rsw increases. This is due to the 
wind-terminating shock becoming highly distorted so the 
time for the wind to travel from the star to the shock varies 
from pole to equator. Therefore if the wind is varying on a 
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Figure 2. Density profiles through the simulations shown in Fig- 
ure [JJ The solid and dotted lines are through the spherically- 
symmetric wind while the dashed and dash-dotted lines are 
through the asymmetric simulation along the pole and equator 
directions respectively. 



similar timescale the wind-termination shock will respond 
to any change on a different timescales. 

With these distances and wind densities it is possible to 
estimate how long it would take for a GRB jet to traverse the 
distance to the beginning of the stalled- wind at Rsw- Rather 
than perform a detailed afterglow calculation we use a sim- 
ple model of isotropic relativistic ejecta with Ei BO — 10 53 ergs 
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Figure 3. Similar to Figure [T] but 340,000 years later, a few years before core-collapse. The left panel shows the case for a spherically 
symmetric wind while the right panel shows the case when the equator-to-pole wind density ratio is 5. The polar direction is along the 
x-axis in both figures. 
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Figure 4. Density profiles through the simulations shown in Fig- 
ure \3\ The solid and dotted lines are through the spherically- 
symmetric wind simulation while the dashed and dash-dotted 
lines are through the asymmetric simulation along the pole and 
equator directions respectively. 



and Finitiai = 100. As the ejecta expands it sweeps up mate- 
rial decreasing it's velocity. Therefore shorter distances and 
lower wind-densities both shorten the time for the jet to 
reach Rsw- We note that we calculate the observer time 
that includes the effect of time dilation by viewing along 
the expanding ejecta direction of motion. However it does 
not include cosmological time-dilation therefore these times 
are a lower limit and will be longer for higher redshifts. 



Table [T] show that the asymmetries greatly reduce the 
time for the jet to reach Rsw- However a high density ISM 
is still required to reduce the size of the wind-bubbles. The 
time of around an hour is short enough that the early X-ray 
afterglow of GRBs may show some sign of occurring in a 
free-wind region while later emissi on may show a c o nstant 
density profile. This agrees with IPanaitescu et all (|2006l ) 
where for the first few hours of a sample of GRB after- 
glows a free- wind CBE is preferred while after a few days 
a con stant density CBE become more common l|Panaitescul 
2005). Therefore we may see the transition from free-wind 
to stalled- wind regions in some GRB afterglows without any 
bump in the lightcurve when the jet encounters the wind 
termination shock. It is important to note that our time es- 
timate is approximate, if for example the Lorentz factor or 
the kinetic energy of the ejecta is higher then the times can 
be reduced further. Although our results demonstrate that 
even moderate asymmetries in the wind reduce the times by 
an order of magnitude. 



4 DISCUSSION AND CONCLUSIONS. 

While the asymmetry in WR winds of WR151 and other 
WR stars could be due to rotation, duplicity or some other 
factor, it is the asymmetry itself that has a large effect 
on the stellar-wind bubble structure. It may also affect 
the evolution of the star since angular momentum loss will 
depend on how much m ass is lost at equatorial latitudes 
<|Mevnet fc Maederl l2007h . It will effect the evolution of a 
SN or GRB if it occurs within such a distorted wind-bubble. 
Spherically symmetric SN ejecta will be forced into an oblate 
form. For a GRB if we assume that the observed asymme- 
try is caused by rotation and that the direction of the low- 
est wind density is the rotation axis then the GRB jet will 
propagate in that direction. Therefore the more asymmetric 
the WR wind the more quickly it will encounter the wind- 
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Table 1. Position of the wind-termination shock from the WR 
stars with different initial ISM densities and different equator- 
to-pole density ratios during the WR phase. Also listed is the 
wind density in the polar direction for the different density ratios, 
A» = (A//10~ 5 A/ Q yr^V^wind/lOOOkms- 1 ). Also show is the 
time taken for ejecta to reach the wind- termination shock -Rgw for 
a simple isotropic model with E\ so = 10 53 ergs and r; n ; t ; a i = 100. 
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In conclusion some WR stars are observed to have asym- 
metric winds and this has a large effect on the position of 
the wind-termination shock. The reduction of the wind den- 
sity in one direction can lead to the wind-termination shock 
moving closer to the star. If a GRB jet was to travel in this 
direction then with the lower wind density and shorter dis- 
tance, a constant density CBE is more likely to be inferred 
from the observations of the afterglow. 
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termination shock taking only a few hours. Furthermore be- 
cause the wind density is also reduced the jet will also be 
decelerated less as it will sweep up less material. 

These factors together indicate that the afterglow jet 
is more likely to reach the constant density environment 
quickly so the free-wind environment does not affect the 
afterglow. This is suggested by the times listed in Table 
[T] Combined with the oth er possible mechanisms discussed 
by Ivan Marie et al.l (|2006T ) to reduce the wind-termination 
shock distance we should not be surprised that many GRB 
afterglows apparently occur in constant density environ- 
ments rather than free-wind environments. Also there is 
some evidence that stellar wind mass-loss rates are lower 
than currently used in stellar mode ls with stars losing mos t 
of their mass in explosive outbursts l|Smith fe Owockil 2006). 
If mass-loss rates and therefore wind densities are reduced 
again a constant density environment becomes more likely. 

We must also consider that GRB progenitors may re- 
quire very rapid rotation. If this is the cas e then the mod- 
els of llgnace. Cassinelli fe Biorkmanl (1 19961) ma y not a pply 
as pointed out by Qwocki. Cranmer fe Blondinl (|l994T ) and 
lOwocki. Cranmer fc Gavlev 1 19961 ). If this is the case the 
reverse of our argument will be true and the wind termina- 
tion shock distance along the polar direction will be much 
greater and the wind much denser. It is difficult to know ex- 
actly which model is correct, not because the models of the 
wind acceleration are limited but rather because the stel- 
lar evolution models are limited by being only ID models. 
In the case of rapid rotation, the stellar structure will be- 
come distorted and ID models give some indication of what 
might occur but they cannot fully model this inherently 3D 
problem. 

By using the observed results on WR winds of what 
we assume is rotation (or duplicity) then we bypass this 
uncertainty. One remaining problem is that the WR stars 
observed to have asymmetries to date are solar metallicity 
WR stars in the Milky Way and it is though that GRB pro- 
genitors are lower metallicity WR stars. Therefore a study 
into the polarisation and asymmetry of lower metallicity WR 
stars in the LMC and SMC would be of great interest (e.g. 
Vink, in prep.). 
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